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ABSTRACT 

_ We present B and R band spectroastrometry of a sample of 45 Herbig Ae/Be 

stars in order to study their binary properties. All but one of the targets known to be 
binary systems with a separation of ~ 0.1 — 2.0 arcsec are detected by a distinctive 
spectroastrometric signature. Some objects in the sample exhibit spectroastrometric 
Q ■ features that do not appear attributable to a binary system. We find that these may 

be due to light reflected from dusty halos or material entrained in winds. We present 
C/3 , 8 new binary detections and 4 detections of an unknown component in previously 

discovered binary systems. The data confirm previous reports that Herbig Ae/Be 
stars have a high binary fraction, 74 ± 6 per cent in the sample presented here. We 
. use a spectroastrometric deconvolution technique to separate the spatially unresolved 

^ ' binary spectra into the individual constituent spectra. The separated spectra allow 

I us to ascertain the spectral type of the individual binary components, which in turn 

I allows the mass ratio of these systems to be determined. In addition, we appraise 

' the method used and the effects of contaminant sources of flux. Wc find that the 

distribution of system mass ratios is inconsistent with random pairing from the Initial 
Mass Function, and that this appears robust despite a detection bias. Instead, the 
mass ratio distribution is broadly consistent with the scenario of binary formation via 
0^ . disk fragmentation. 
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Key words: binaries: general - stars: emission-line - stars: pre-main-sequence - bi- 
naries {including multiple): close - techniques: spectroscopic 



1 INTRODUCTION derstood (|Bouvier et all 120071 ). In the case of more massive 

^ , ,. r 1 r ■ n, ,.r stuTs the situatioH is much less clear. Such stars do not ex- 

Our understandme ot the tormation and early evolution ot . j_- n • -i i nivir. i i i 

X, penence an optically visible FMb phase, evolve on a much 
massive stars (M* > 8Mn ) is much less complete than in the . , . . , , • i i i i • 

^ ~ r, r more rapid timescale, and are considerably more luminous 

case or low mass stars. The scenario ol low mass star torma- 
tion has been relatively well studied, and a broadly consis- 



than low mass stars. Early studies on the effects of radia- 
tion pressure and the considerable ionising output of massive 
tent observational and theoretical picture has now emerged. , , , i i- j_i ^ ■ j_ r 

, , , , young stars prompted speculation that massive star torma- 
The various phases oi low mass star tormation include: cloud ,. . ,. , . . n 

tion might proc eed m a different manner to t hat ot low mass 
coUapse, proto-steUar creation and a subsequent contrac- ^^^^^ (Larson fc Starrfieldim E^\mi)- For example, 
tion of Pre Main Sequence (PMS) obiects towards the Zero -j, i i_ ^ i j_i ^ j_i ^ ■ j_ i 

, , . „ ,^ . , rrnn . , , m . ^ tias beeu suggested that the most massive stars torm via 

Age Mam Sequence (ZAMS . This later stage, the T Taun , , . ITTT, — — I 

, , , r 11 stellar mergers or competitive accretion (IBallv fc Zmneckeii 

phase, is easy to observe and therefore relatively well un- [200d) 



Based on observations made with the William Herschel Tele- 



However, recent work, on both the observational and 
theoretical front, suggests that massive star formation may 

, , , ^ , , , . , , r not be dissimilar to low mass star formation. As an exam- 
scope and the Isaac Newton telescope operated on the island ot i r i FT^ \ Tl \ Is 

La Palma by the Isaac Newton Group in the Spanish Observato- observational resuffs, | Patel et al. | report the 

rio del Roque de los Muchachos of the Instituto de Astrofisica de detection of a massive disk around a 15Mg protostar, in- 

Canarias. dicating that massive stars may form via monolithic ac- 

t E-mail: pyhewOleeds. ac.uk cretion. On the theoretical front, recent work indicates ac- 
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cretion onto a massive protostar is not impeded by radia- 
tion pressure (lYorke fc Sonnhalteij[20o3 : [Turner et allbOOTl : 
iKrumholz et all 120091 )" However, while significant progress 
has been made, there remain many unaddressed questions 
related to the formation and evolution of massive stars 
ijZinnecker fc Yorkel2007l 'l. As observations of massive young 
stars are challenging, the full extent of the differences and 
similarities between low and high mass star formation are 
still unknown. 

Between the two extremes of mass lie the Herbig Ae/Be 
(HAe/Be) stars (|Herbiel Il960l ). These stars represent the 
most massive of objects to experience an optically visible 
PMS evolutionary phase. Therefore, HAe/Be stars offer an 
opportunity to study the early evolution of stars more mas- 
sive than the sun. Spectropolarimetry indicates that Herbig 
Ae stars may undergo a PMS phase similar to that of the T 
Tauri stars, while Herbig Be stars may evolve via disk ac- 
cretion, rath e r than magnetospher ic accretion l|Vink et al.l 
|2002| . l2005al : iMottram et"al] l2007h . Therefore, it appears 
that a transition in formatio n mechanisms occurs across 
the HAe/Be mass boundary l|Mottram et al.|[2007h . How- 
ever, the critical mass has not yet been established. 

To examine the similarities and differences between low 
mass T Tauri stars, HAe/Be stars and the optically in- 
visible Massive Young Stellar Objects (MYSOs), study of 
the circumstellar environment at small angular scales is re- 
quired. This is not trivial, requir ing observations w ith high 
angu l ar resolution I Mannings fc Sarg ent 1997; Fuc nte et al.l 
I2OO6I : lGradvll2007l : iKTaus et al.1 120081 ). Despite the progress 
in the field, a full understanding of HAe/Be stars is ham- 
pered by the s mall s ample sizes involved. By way of contrast, 
iBaines et al.l ('2006) utilised spectroastrometry to study a 
large sample of HAe/Be stars with milli-arcsecon d (mas) 
precision. Despite this resolution iBaines et al.l l|2006l ) did not 
detect any accretion disks around HAe/Be stars. However, 
they did find that the majority, 68 ± 11 per cent, of HAe/Be 
stars reside in relatively wide (probably a few-hundred au, 
see Section [5. 3p binary systems. 

The binary fraction reported bv lBaineT et all l|2006l l IS 
greater than that of T Tauri stars at similarly wide separa- 
tions, which in turn is greater than t hat of Main Sequence 
G-dwarfs at the same separations (IDuauennoY fc Mavorl 
I1991I : iReipurth fc Zinneckerlll993l : iGhez et al.ll 19931 ). Indeed, 
this high b inary fraction is appr oaching that of more mas- 
sive stars (jPreibisch et al.lll999l ). However, little is known 
about the properties of such binary systems. The proper- 
ties of the binary components and configurations of such 
systems are of interest as they can constrain the binary for- 
mation mechanism. The se minal study to date is that by 
iBouvier fc CorporonI (|200lh . who used Adaptive Optics as- 
sisted observations to construct Spectral Energy Distribu- 
tions (SEDs) for each component in a number of HAe/Be 
binary systems. The drawback of SED fitting is that PMS 
stars, as young stars, are inevitably associated with dusty, 
obscured environments. Therefore, the brightness ratio of a 
binary determined by SED fitting can occasionally be am- 
biguous. However, very few HAe/Be binary systems have 
been studied with spatially resolved spectroscopy, and thus 
far such st udies have been conducted with seeing limited 
resolution (jCarmona et al.ll2007l : iHubrig et al.ll2007l ). 

The position angles of HAe/Be binary systems seem 



to be preferentially aligned wit h the spectropolari metrically 
detected circumprimary disks l|Baines et al.fl2006l ). This al- 
ready places constraints on the formation modes of these 
stars, in that it seems the systems formed via fragmentation 
of a molecular core or d isk. This had already been suggested 
for lo wer mass binaries l|Wolf et al. l l200ll : lKroupa fc Burkerl 
l200j ). but little is known about the formation mechanisms of 
more massive stars. This pa per describes a spec troastromet- 
ric follow-up of the work of IBaines et al.1 (POOQ ) with dedi- 
cated observations to study both components of binary sys- 
tems. The objective is to determine the properties of these 
binary systems and thus place stronger, more quantitative, 
constraints on the formation of stars of intermediate mass. 
We do this by determining the mass ratio of these binary 
systems. This is done using a spectroastrometric technique 
to disentangle the constituent spectra of unresolved binary 
systems, allowing the spectral type, and hence mass, of each 
component to be determined. Spectroastrometry itself is a 
relatively simple technique that extracts the spatial informa- 
tion present in conventional longslit spectra. Crucially, spec- 
troastrometry can probe changes in flux distributions with 
a typical precision of a mas or less (jBailoy 1998a), which is 
required to study unresolved binary systems. Typically the 
minimum separation probed is of the order 100 mas, as the 
signature of a binary system is dependant upon the system 
brightness ratio and separation. 

This paper is structured as follows: in Section [2] we 
present our sample selection, observation method and data 
reduction procedures. In Section[3]we discuss the spectroas- 
trometric signatures observed. In Section |4] we present the 
method of splitting unresolved binary spectra and in Section 
14.11 we review the results of separating binary spectra into 
their constituent spectra. In Section[5]we discuss our results. 
Finally, we conclude this paper in Section |S] by summarising 
the salient points raised. 



2 OBSERVATIONS AND DATA REDUCTION 
2.1 Observations 

The data presented consist of long-slit spectra in the B band 
(4200 - 5OOOA) and/or the R band (6200 - 7000A) of 45 
HAe/Be stars, and 2 emission line objects which are possi- 
bl e HAe/Be s t ars. The objects were chos en from the catalogs 
of 'T he et al] l|l994h . lVieira et al.1 (|2003l ) fc lHernandez et al] 
(20041 ) ■ and were selected to be reasonably bright ( F ^ 
12-13) . Some objects previously observed by IBaines et al.l 
(!2006|) were observed to provide a consistency check on the 
spectroastrometric signatures. Given the small population 
of HAe/Be stars, the objects observed constitute a repre- 
sentative sample of HAe/Be stars, albeit brightness limited. 

The data were obtained using the 4.2m William Her- 
schel Telescope (WHT) and the 2.5m Isaac Newton Tele- 
scope (INT). At the WHT, data were obtained on the 
6th & 7th of October 2006, using the Intermediate Disper- 
sion Spectrograph and Imaging System (ISIS) spectrograph. 
Spectra of 20 objects were taken simultaneously in the B and 
R bands using the dichroic slide of ISIS. In most cases a slit 
5 arcsec wide was used to ensure all the light from a given 
binary system entered the slit, even in poor seeing. This al- 
lows us to study the individual binary components, unlike 
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iBaines et al] (|2006l '). who used a slit of 1 arcsec. The R1200B 
and R1200R gratings were used and the resulting spectral 
resolving power was found to be ~ 3500, corresponding to 
85kms~^. The angular pixel size was 0.20 and 0.22 arcsec 
in the B band and R band respectively, which means that 
the spatial profile of the longslit spectra was well sampled 
(average FWHM 1.9 arcsec). At the INT data were obtained 
using the 235mm camera and the Intermediate Dispersion 
Spectrograph (IDS). Observing was conducted from the 27th 
of December 2008 to the 3rd of January 2009. The spectra 
of 32 objects were obtained, despite adverse weather condi- 
tions preventing observing for the better part of three nights. 
As at the WHT the slit width was generally 5 arcsec. The 
R1200R and R1200B gratings were used and the resulting 
spectral resolution was found to be ~ 3800, or 80kms~^. 
The angular size of the pixels was 0.4 arcsec, which fully 
sampled the average spatial profile of the spectra (1.8 arc- 
sec). 

Multiple spectra were taken at four position angles (PA) 
on the sky. The PAs selected always comprised of two per- 
pendicular sets of two anti-parallel angles, e.g. 0° , 90° , 180° 
and 270° . Dispersion calibration arcs were made using CuNe 
and CuAr lamps. Table [T] presents a summary of the obser- 
vations. 



2.2 Data reduction 

Data reduction was conducted using the Image Reduction 
and Analysis Facility (IRAF j3 and routines written in Inter- 
active Data Language (IDL). Initial data reduction consisted 
of bias subtraction and flat field division. The total inten- 
sity spectra were then extracted from the corrected data in a 
standard fashion. Wavelength calibration was conducted us- 
ing the arc spectra, and the wavelength calibration solution 
had a precision of the order < O.lA. 

Spectroastrometry was performed by fitting Gaussian 
functions to the spatial profile of the long-slit spectra at each 
dispersion pixel. This resulted in a positional spectrum, the 
centroid of the Gaussian as a function of wavelength, and 
a Full Width at Half Maximum (FWHM) spectrum, the 
FWHM as a function of wavelength. Spot checks were used 
to ensure that a Gaussian was an accurate representation of 
the data. The continuum position exhibited a general trend 
across the CCD chip: of the order of 10 pixels in the case of 
the ISIS data and 2 pixels in the IDS data. This was removed 
by fitting a low order polynomial (4**^ or 5**^ order) to the 
continuum regions of the spectrum. 

All intensity, positional and FWHM spectra at a given 
PA were combined to make an average spectrum for each 
PA. A correction for slight changes in the dispersion across 
PAs was determined by cross-correlating average intensity 
spectra obtained at different PAs. The correction was then 
applied to the average intensity, positional and FWHM spec- 
tra. The average positional spectra for anti-parallel PAs were 
then combined to form the average, perpendicular, position 

^ IRAF is distributed by the National Optical Astronomy Obser- 
vatories, which are operated by the Association of Universities for 
Research in Astronomy, Inc., und er cooperat ive agreement with 
the National Science Foundation llTodvlll993l) . 



spectra, for example: (0° - 180° )/2 and (90° - 270°)/2. This 
procedure eliminates instrumental artifacts as real signa- 
tures rotate by 180° when viewed at the anti-parallel PA, 
while artifacts remain at a constant orientation. In addi- 
tion, all positional spectra were visually inspected for arti- 
facts not fully removed by this procedure. As with the posi- 
tional spectra, the FWHM spectra at anti-parallel PAs were 
also combined to make two averaged, perpendicular spectra. 
While FWHM features do not rotate across different PAs, 
the features observed at anti-parallel PAs were used to ex- 
clude artifacts via a visual comparison. All conditions being 
constant, a real FWHM signature should not change from 
one PA to the opposite angle at +180°. 



3 SPECTROASTROMETRIC SIGNATURES 

3.1 Binary spectroastrometric signatures over Hi 
lines 

An unresolved binary system, in which each component has 
a unique spectrum, displays a clear signature in the be- 
haviour of the spectral photocentre. As the spatial profile 
is the sum of the two stars convolved with the seeing, the 
peak is not located at the position of either star, but some- 
where between the components. The exact location of the 
photo-centre depends on the intensity ratio and the separa- 
tion of the two components. Over spectral lines the binary 
flux ratio changes from its continuum value, which results 
in the peak position shifting towards the dominant compo- 
nent. Therefore, unresolved binary systems are revealed by 
a displacement in the positional spectrum over spectral line. 
In addition, an unresolved binary system is also revealed by 
a change in the FWHM over lines in the spectrum. Again, 
this is because the spatial profile of an unresolved binary sys- 
tem is dependent upon the binary fiux ratio, which changes 
from its continuum value across certain lines. As the error 
in the centre of the Gaussian profile is governed by photon 
statistics, changes of mas scales can be traced. This allows 
binary systems with separations as small as ~ 0.1 arcsec and 
differences in brigh tness up to 5 magnitudes to be studied 
(|Baines et al.ll2006D . 

To illustrate the detection of a binary system the obser- 
vations of GU CMa are presented in some detail. GU CMa 
is known to be a Herbig Be binary system with a separation 
of 0.65 arcsec, a PA of 195°, a brightness difference be- 
tween components of 0.7-1.0 magnitudes in t he optical band 
and a primary with a spe c tral type of Bl ('Fu et al.l ll997l : 
iFabricius fc Makarovllioool : iBouvier fc Corp oron 200 3). 

GU CMa presents a very clear binary signature in the 
spectroastrometric observations (Fig.[l|. Across the Hi lines 
the photo-centre of the spectrum clearly shifts towards the 
North-East. This demonstrates that the primary, the com- 
ponent brightest in the continuum, dominates the emission 
spectrum. It also indicates that the secondary, the compo- 
nent least bright in the continuum, has the larger absorption 
profile of H7. As the photo-centre shifts to the North-East 
the FWHM of the spectrum is seen to decrease. This also 
indicates that the primary dominates the spectrum at these 
particular wavelengths. The photo-centre is also observed to 
shift to the North-East across the He I lines. This again in- 
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Table 1. Log of the observations, column 1 lists the objects observed, column 2 denotes the spectral type of the objects, column 3 lists 
the V band magnitudes of the sample, and column 4 designates which telescope the object in question was observed with. Columns 5 
and 6 list the average seeing conditions, columns 7 and 8 list the total exposure times and column 9 denotes the slit width used. Column 
10 lists the total Signal to Noise Ratios, and finally, column 11 presents the date(s) each object was observed. Information on the objects 
is taken from SIMBAD (simbad.u-strasbg.fr) unless otherwise stated. 



Object 


Spec type 


V 


Telescope FWHM'' EWnM*" 


^bluc 


^rcd 


Slit 


SNR 


Date 










{") 


{") 


(s) 


(s) 


( ) 






V -A. V^dyO 


AOe 


1 1 


WHT 


1.3 


1.2 


Asnn 

4:OUU 


^ouu 


^ n 
o.u 


UUUB lO / ur 


07/1 n /9nnfi 
u ( / -LU/ zuuu 


V -A. '^da 


/\ue 




INT 


1.1 


1.7 


9snn 
zouu 


ODUU 


Z.OB) 


O / UB 1 O ' UR 


9!R/i9/9nn!R "^i /1 9 /9nns 

ZO/ IZ/ZUUOiOl/ iz/zuuo 


V Ot}'± v^jdb 


Be 


J-U-U 


INT 


1.3 


- 


'?9nn 




^ n 
o.u 


fil n 

UlU 


ni /ni /9nnQ 
U-L/U-L/ zuuy 


\T^ 1 Tail 
V lloO -LciU 


A 1 


10 7 
iU. 1 


INT 


1.7 


- 


OZUU 




o.u 


^OU 


fiQ /ni /9nnQ 
uo/ ui / zuuy 




A '^ 


1 n 

iU.O 


INT 


1.2 


1.4 


9nnn 
zuuu 


9Ann 
z^uu 


z.OB) o.Ur, 


i lUB jO'^UR 


9!R/i9/9nn!R "^i /1 9 /9nns 

ZO/ ±Z/ZUUO,Ol/ IZ/ZUUO 




rY\J V pe 


7 1 


WHT 


1.9 


1.9 


ooU 


ozu 


o.u 


1 1 n^i fii^nx-. 

i ±Ub )DoUr 


c\(\ /in /9nn(^ 

UD/ ±U/ ZUUD 


iVi VV K-f 4:OU 


A3pshe 


7 7 


WHT 


2.0 


2.1 


yuu 


u^u 


^ n 
o.u 


1 1 nnn Rnn^ 

IIUUB )0UUR 


r\(\ /I n /9nnfi 
uu/ J-U/ zuuu 


UX Ori 


A3e 


y .u 


WHT 


2.4 


2.4 


ouuu 


ouuu 


^ n 
o.u 


1 9nnx-. QAnx-, 

1 ZUUR ) y^^R 


c\(\ /I n /9nnfi 
uu/ J-U/ zuuu 


\ L\J Wi 1 




12 1 


INT 


2.1 


- 


Asnn 

4:OUU 




^ n 
o.u 


1 

lOU 


09 /ni /9nnQ 
uz/ ui/ zuuy 


V -LOUU \Jl 1 


AOe 


y .o 


INT 


1.3 


- 


9Ann 




n 

o.u 


^7n 

O ( u 


Qi /1 9 /9nnR 

O-L/ -LZ/ ZUUO 


\T'iA(\ Ori 
V O^iJ \Jl\ 


A ^TTT 


1 n 1 


INT 


1.5 


- 


ODUU 




o.u 


9nn 
zuu 


ni /ni /9nnQ 
ui/ Ui/ zuuy 




A ^ 


S 1 

O. -L 


WHT 


1.9 


1.7 


9nfin 
zuuu 


1 A7n 


^ n 
o.u 


AC\-r^ Qnnx-, 
"y^UBiyiJUR 


07/1 n /9nnf^ 

U ( / -LU/ ZUUU 


V OoU Wi 1 


An 


10 7 
lU. / 


INT 


1.5 


1.6 


ouuu 


9Q4n 
zy4:U 


Q Or-. ^ On 

O.UB ) O.UR 


1 nnn 9nnTi 

lUUB j^'JDR 


9)R /1 9 /9nn)R "^i /i 9 /9nn^i 

ZO/ -LZ / ZUUO, O-L/ -LZ / ZUUO 


iVi VV Kj I oo 


A3e 


o.o 


WHT 


1.4 


1.3 


-LUOU 


yuu 


^ n 
o.u 


OUb i^'J'JR, 


n7 /I n /9nnf^ 
u / / J-U/ zuuu 


HK Ori 


.A4pev 


1 1 Q 
1 ±.y 


INT 


2.4 


— 


Asnn 

lOUU 




^ n 
o.u 


9nn 

zuu 


09 /ni /9nnQ 
uz / u J- / zuuy 


n.±y ^'-t'-tyjyj'-t 


A3 


9.4 


WHT 


1.7 


1.6 


3180 


3660 


5.0 


inn^ 79nTi 

J-UUB ) ' -^UR 


07/1 n/90nR 
u / / J-U/ zuuu 


V J-Z ( ± WI 1 


A 


inn 


INT 


1.6 


- 


z^uu 




^ n 
o.u 


410 


01 /Ol /90nQ 
ui / u J- / zuuy 


T Ori 




y . o 


INT 


1.9 


— 


'?9nn 
ozuu 




^ n 
o.u 


•^nn 

OUU 


ar> /ni /9nnQ 

uo 1 i- / zuuy 


■\/"r:Qf; Ori 


A2V 


y .o 


INT 


3.3 


— 


9Q4n 

zyiu 




^ n 
o.u 


uou 


09 /ni /9nnQ 
uz / u J- / zuuy 


nu o 1 oo I 


AOe. 


o.o 


INT 


1.4 


1.4 


9nfin 
zuuu 


1 A7n 


O.UBi O.Ur 


ZUUB jOOUR 


9>^ /1 9 /9nn)R "^i /1 9 /9nn^i 

ZO/ J-^/ ZUUO, O X/ J-Z/ ZUUO 


\r-\ 7QQ Ovi 
V 1 1 oo KJTl 


RQA/o 

-Dy V c 


Q Q 

y .y 


INT 


1.7 


- 


lOOU 




o.u 




01 /ni /9nno 
ui/ Ui/ zuuy 




13 yi V 


1 n 7 

-LU. / 


INT 


1.8 


— 


9snn 
zouu 




^ n 
o.u 


1 nn 

lUU 


0*^ /ni /9nnQ 
uo/ Ui/ zuuy 


RR Tan 


A9TT TTTt. 


1 Q 

lu.y 


INT 


1.7 


— 


9snn 
zouu 




^ n 
o.u 


9'^n 

zou 


0*^ /ni /9nnQ 
uo / u i / zuuy 


"\/"QC:n Ori 
V OUU WI 1 


AOe 


^(^ A( fi^ 


INT 


1.9 


— 


Asnn 

4:OUU 




^ n 
o.u 


1 'in 

lOU 


0"? /ni /9nnQ 
uo/ Ui / zuuy 


M wo 1 9n 

iVl VV _LZU 


An 


7 Q 

/ . y 


WHT 


2.1 


1.9 


4sn 

4:OU 


^ou 


^ n 
o.u 


louuB )UyuR 


c\f\ /1 /9nnfi 
uu/ iu/ zuuu 


M wo T 9n 

iVi VV Kj -LZU 


An 


7 Q 

1 . y 


INT 


1.4 


1.6 


z^uu 


1 9^n 


O.UB) o.Ur 


IZUUR i^'J^R, 


9S /1 9 /9nn>^ "^1 /1 9 /9nn^i 

ZO/ i-^/ ZUUO, O i/ iZ / ZUUO 


MWO 7Qn 


Be 


12 


INT 


3.1 


— 


4:UOU 




^ n 
o.u 


9nn 

zuu 


09 /ni /9nnQ 
uz / Ui / zuuy 


MWP 1 '^7 

iV J. w \^ ±o 1 


Be 


11.2 


INT 


— 


2.1 




4560 


5.0 


200 


9)^ /I 9 /90n^^ 

ZO/ iZ/ ZUUO 




Bpshe 


S n 
o.u 


WHT 


2.0 


1.9 


OUU 


240 


^ n 
o.u 


yuuB jOOUR 


ofi /I n /9onfi 
uu/ iu/ zuuu 


T,VHrv 91 ^ 

i_J IV± Z( J. o 


B7.5e 


10.6 


INT 


1.5 


2.3 


3600 


3600 


O.UB ) ^-^R 


ouuB )OUUr 


97/19/900^^ "^l /19/90n^^ 

Z / / iZ / ZUUO, O ±/ ±Z / ZUUO 


MWO 1 A 7 

iVi VV -L'i 1 


B6pe 


o.o 


WHT 


1.8 


1.5 


"^nnn 
ouuu 


1 7nn 

1 1 uu 


^ n 
o.u 


1 9nnx-. f^nnxi 

IZUUB WJUR 


07/1 n /90nfi 
u / / iu/ zuuu 


MWP 14-7 

iVl VV V_' -Lrt ( 


B6pe 


8.8 


INT 


1.3 


— 


2400 




5.0 


620 


ni /Ol /9onQ 

ui/ui/ zuuy 


LX, iViUll 


RO 

-DU 


10 4 


INT 


4.2 


2.5 


ouuu 


ouuu 


^ n 
o.u 


1 nnn '7AC\r, 

lUUB j^^DR 


9!>i /I 9 /90n)R n9 /Ol /90nQ 

ZO/ iz / zuuo,uz / Ui / zuuy 


■\/"c;Qn A/Tnn 

V oyu iviuii 


B8pe 


1 9 Q 
iz.y 
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1.5 




9fi7n 
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^ n 
o.u 


9nn 

zuu 


01 /Ol /900Q 
ui/ui/ zuuy 


V742 Mon 


B2Ve 


6.9 


INT 


1.4 


2.8 


1740 


2535 


5.0 


400b,800r 


30/12/2008,31/12/2008 


OY Gem 


Bp[e] 


11.1 


INT 


1.8 




2880 




5. 


100 


03/01/2009 


GU CMa 


B2Vne 


6.6 


WHT 


2.5 


2.4 


360 


360 


5.0 


1500b, 900r 


06/10/2006 


GU CMa 


B2Vne 


6.6 


INT 


1.8 




720 




5.0 


1400 


03/01/2009 


MWC 166 


BOIVe 


7.0 


WHT 


2.5 


2.3 


210 


120 


5.0 


1200b, 800r 


06/10/2006 


HD 76868 


B5 


8.0 


INT 


1.5 


2.4 


4830 


2100 


5.0 


100b, IOOr 


30/12/2008,01/01/2009 


HD 81357 


B8 


8.4 


INT 


1.7 




4800 




5.0 


100 


03/01/2009 


MWC 297 


Be 


12.3 


WHT 


1.4 


1.2 


4100 


3120 


5.0 


100b,500r 


07/10/2006 


HD 179218 


B9e 


7.2 


WHT 


2.6 


2.4 


2100 


1200 


1.0/1.5 


2300b, 940r 


06/10/2006 


HD 190073 


A2IVpe 


7.8 


WHT 


1.5 


1.2 


540 


360 


5.0 


600b,800r 


07/10/2006 


BD +40 4124 


B2 


10.7 


WHT 


2.1 


1.9 


600 


660 


4.0 


500b,370r 


07/10/2006 


MWC 361 


B2Ve 


7.4 


WHT 


1.7 


1.7 


1350 


960 


2.5/4.0 


1400b, 1400r 


06/10/2006 


SV Cep 


Ae 


10.1(B) 


INT 


1.6 




3000 




5.0 


600 


02/01/2009 


MWC 655 


BllVnep 


9.2 


INT 


1.7 




2400 




5.0 


400 


03/01/2009 


11 Cep 


B2IV/Ve 


9.3 


WHT 


1.4 


1.2 


3500 


3000 


5.0 


800b,500r 


07/10/2006 


BHJ 71 


B4e 


10.9 


WHT 


1.8 


1.8 


1200 


1080 


4.0 


500b,340r 


06/10/2006 


BHJ 71 


B4e 


10.9 


INT 


1.7 




4200 




5.0 


500 


01/01/2009 


MWC 1080 


BO 


11.6 


WHT 


2.0 


2.0 


3300 


4170 


5.0 


200b,500r 


06/10/2006 



^ Average seeing in the blue spectral region, approximated by the average of the individual median FWHM, where necessary averaged 
over multiple slit widths. 

^ Average seeing in the red spectral region, approximated by the average of the individual median FWHM, where necessary averaged 
over multiple slit widths. 
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dicates that it is the primary that dominates the binary flux 
over this hne. 

AU the spectroastrometric excursions across the H I 
hues produce a PA consistent to within ~ 1°. The agree- 
ment between the spectroastrometric displacements across 
different hues is an important consistency check. In addition, 
the close agreement between the data gathered at different 
telescopes provides compelling evidence that the signatures 
observed are real and not contaminated by instrumental af- 
fects. The difference in the FWHM changes can be explained 
by the difference in the seeing between observat ions. Finally, 
that our observations concur with the results o fH aines et al.l 
further proves that these signatures are real and no 
instrumental effect. As demonstrated by the spectroastro- 
metric signature associated with the He I line such features 
do not solely occur across lines with an emission component. 
This is important to note as it means the spectra splitting 
method does not require emission lines to separate spectra. 

In the case of some stars signiflcant FWHM changes 
are observed which ar e not accompanied by a change in the 
spectral photo-centre. iBaines et al] (|2006l ) regard such a sig- 
na ture as a p ossible binary detection. This is substantiated 
as iBaines et al. (^006) demonstrate that the spectroastro- 
metric signature of a binary system with a separation of 
greater than half the slit width exhibits larger FWHM than 
positional features. However, in the data presented here, this 
scenario is unlikely. As a wide slit was used, a binary with a 
separation of half the slit width would be resolved, even in 
seeing conditions of 2 arcsec. As a resolved system does not 
exhibit a spectroastrometric signature we suspect there may 
be an alternative explanation to the FWHM features not 
accompanied by positional features. We note that such fea- 
tures are not instrumental as some stars exhibit no change in 
FWHM over spectral lines. In many cases the large FWHM 
features occur over absorption features in the emission pro- 
files. This suggests that these features may trace an extended 
structure which scatters the line profile, rather than being 
an intrinsic source. If the scattering media were close to 
being symmetrically distributed around the central star it 
could generate a large FWHM increase while not resulting 
in a positional si gnature. Such source s of flux could be a 
dis k/stellar wind (lAzevedo et al] 2007^■ th e halos reported 
bv lLeinert et all l|2001^ and lMonnier et~al] (|2006l '). or nebu- 
losity. This topic will be returned to in Section [5. II 

The observational results naturally fall into three cat- 
egories: clear binary signatures, possible binary/other sig- 
natures and null detections. We present a summary of the 
detections in our results in Table (2] in which we separate 
known binary systems and new spectroastrometric detec- 
tions. It is important to note that it is not only binaries 
that are detected by spectroastrometry, optical outflo ws and 
disks can also result in a s pectroastrometric feature l|Bailevl 
ll998bl : lTakami et al.|[200ll ). However, there are no disk sig- 
natures, and only a few detections of outflows, in the data 
presented here. 

There are 29 stars in our sample that are referred to 
in the literature as being part of a binary system. However, 
we exclude AB Aur and HD 244604 as the binary nature of 
these objects is open to question, section 15.11 explains why 
this is the case. We detect 20 of the 27 previously known bi- 
nary systems. As six of the seven undetected systems have 



separations greater than 2 arcsec and/or brightness differ- 
ences as great as 8 magnitudes, we detect all but one, UX 
Ori, of the binary systems that we would expect to detect. 
Given that the majority of known binary systems are de- 
tected, objects which are not known to be part of a binary 
system but exhibit similar spectroastrometric signatures to 
the known binary systems are classified as new binary de- 
tections. We detect 8 new binary systems. The raw binary 
fraction of the sample is 0.60. Including the non detections 
of known binary systems the binary fraction of the sample 
is 0.74. While these figures are high for a limit ed separation 
range, they are consisten t with previous work l|Pirzkal et al.l 
1199^ IBaines et al.]|2006l ) . 

As expected, we do not detect the binaries with separa- 
tions greater than ~ 2.0 arcsec and differences in brightness 
greater than 5 magnitudes, e.g. HD 179218 and MWC 297. 
We note that the wide companions are not detected in the 
longslit spectra as distinct sources. VX Gas, T Ori, LkHa 
215 and II Gep are all known to be binary systems and all 
display a binary signature in their spectroastrometric sig- 
natures. Therefore, these stars are classified above as detec- 
tions of known binary systems. However, these systems are 
wide binaries with separations greater than 5 arcsec. These 
companions are clearly resolvable, and thus the spectroas- 
trometric signatures we observe can not be due to the pre- 
viously reported companion. Therefore, we suggest we have 
detected previously unknown companions to VX Gas, T Ori, 
LkHa 215 and II Gep. 

Table IXT1 summarises the spectroastrometric signatures 
over Hq and H/3. We note that similar behaviour was ob- 
served across H7 and other lines, but to keep this paper 
concise we only present the Ha and H/9 signatures. Further- 
more, the spectroastrometric signatures of the entire sample 
over either the Ha or H/3 lines are presented in Appendix iBl 
Spectral variability is a common behavi our of H Ae/Be stars 
(|Rodgers et all I2OO2I : iMora et al.ll2004h . In the case of the 
few objects observed twice, some line profile variations are 
seen. However, the spectroastrometric signatures of objects 
observed twice are generally consistent, e.g. the example of 
GU GMa (Fig[T]|. In addition, in the case of objects com- 
mon to the this sample and that of IBaines" et al.l (20061) ■ the 
spectroastrometric signatures presented here are consistent 
with the previous results. Therefore, we conclude that spec- 
tral variability, on timescales of years, does not effect the 
spectroastrometric signatures observed. No line profile vari- 
ability on timescales of minutes is observed. 

To summarise, spectroastrometry is a powerful tool 
with which to study binary systems, as GU GMa demon- 
strates. Not only do we clearly detect a 0.6 arcsec binary in 
seeing as large as 2.5 arcsec, we also trace the PA of the sys- 
tem with a precision of 1° or less. Indeed, spectroastrometry 
detects all but one of the known binary systems with separa- 
tions less than ~ 2 arcsec and differences in brightness of less 
than 5 magnitudes. In addition, the PAs of these systems are 
all traced with a precision of the order of 1°, and are gen- 
erally consistent with literature values to within ~ 5°. Most 
importantly, the spectroastrometric displacements contain 
information as to which component of the binary system 
dominates the flux over certain spectral features. This in- 
formation can be used to separate the constituent spectra. 
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Figure 1. The spectroastromctric signature of GU CMa. In the top panel, presented from left to right: the Ho:, H/9, H7 and He I 
A4471 spectral profiles and associated spectroastrometric signatures. In the spectroastromctric signatures North and East are positive. 
In the lower section of the figure we present the associated XY-plots of the Ha, H/3, H7, and by way of contrast, the He I A 4471 
spectroastrometric displacements. In the XY-plots North is up and East is to the left. The data from the WHT is represented by the 
solid lines, while the data obtained at the INT is represented by the dashed lines. Note the consistency between the two datasets. 
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Table 2. Previously known binary systems and new detections. Column 1 list the objects in question, columns 2 and 3 list the separation 
and PA of the known binaries, taken from the literature. Column 4 contains difference in brightness between the two binary components. 
The amalgamation of data is not complete, if more than one value of binary parameter is available in the literature, only one is presented 
for the sake of clarity. 



Object 


Separation 


PA 


Aflux 




(mas) 


( ) 


(magnitudes) 


Known binaries detected: 






Va Cas 


5340^ 


165.3^ 


/i';4.8* 


V380 Ori 


125 ± 25^ 




K'l 42'^ 


HK Ori 


347.7 ± 2. 5'^ 


41.8 di 0.7^ 


V: 0.87" 


T Ori 


7700 ± 200'^ and spectroscopic^ 


72.6^ 


K:> 4.5*^ 


V586 Ori 


990^ 


3Q 3A 


if ;2.8* 


HD 37357 


186'= 


49 qE 


A" : 1 . 7* 


V1788 Ori 


520^ 


352. 9^ 


7^:3.5* 


HD 245906 


130^ 


77 A 


A';1.5* 


V350 Ori 


290^ 


206.8^ 


A';3.2* 


HD 45677 




150 di 17^ 




LkHa 215 


8500^ 


226.6^ 


A"; 4. 8* 


MWC 147 


150^ 


55 gA 


if :3.8* 


R Mon 


670'^ 


900 


T<'-A 


GU CMa 


654H 


194.5" 


F:0.95±0.02" 


MWC 166 


654H 


297.8" 


y:1.4l" 


BD +40 4124 


720^ 


175.1* 


if;5.4* 


MWC 361 


2250 ± 240' 


164.0 ± 1.0' 


if;4.9' 


SV CepJ 








11 Cep 


6960' 


147.0' 


if:0.0' 


MWC 1080 


760 ± 2<^ 


267.0 ± 1.0<^ 


if;3.25±0.08'' 


Known binaries not detected: 






UX Ori 


22(min)'^ 


257.4 ± 18.4'^ 




MWC 758 


2280^ 


311.3* 


if:8.3* 


V1271 Ori 


8380^ 


294.7* 


if 6. 7* 


V590 Mon 


5007^ 


97.1* 


if6.6* 


MWC 297 


3930 ± 200'^ 


313 ± 2^ 


ii:8.5 ± 0.25!^ 


HD 179218 


2540* 


140.5* 


if:6.6* 


BHJ 71 


6170* 


29.2* 


if8.3* 



New spectroastrometric detections: 

V1366 Ori, HD 35929, RR Tau, MWC 120, V742 Mon, OY Gem, HD 76868 and HD 81357 



References: A) [Thomas et al.l | |2007|| : Bl lSmith et al.l 1120051); O iLeinert et al.l lll997^: D) IShevchenko fc Vitrichenkol lll994l '): 
E) [HartkoDf et al.' ('l99^; F)] Baines et al 

T^m); Cl lWeigelt et al.l ll2002h: H ) 'Fabricius & Makarov" ('2OOOI); 
D IPirzkal et al.. (.1997') : J) Rodgers et al. priv. com. f2008): K'l lBertout et al. (1999); L) Vink ct al. (2005tj). 



3.2 Artifacts 

Several spectroastrometric signatures presented in Table [XT] 
are referred to as artifacts. An artifact is defined as a spec- 
troastrometric signature which does not rotate by 180° when 
viewed at two anti-parallel position angles. Artifacts in spec- 
troastrometric data can arise from a number of sources. In- 
strumental effects include: the misalignment of the disper- 
sion axis with the CCD columns, a change in focus along 
the slit, curvature of the spectrum and any dep arture of 
the CCD array from a regular grid l|Bailevlll998"bl ). As dis- 
cussed in Section[23I such artifacts are readily identified and 
negated. The observation of unresolved lines can also result 
in false signatures (iBa ilcv 1998a). Here, the data were ob- 
tained with a resolution sufficient to resolve most spectral 
lines. However, the narrow absorption troughs in many of 
the double-peaked Hi emission profiles may cause some of 
the artifacts as these features are often barely resolved. 

In addition, iBrannigan et al.] (|2006l ) report an artifact 
that is a consequence of image distortion, regardless of 
whether the spectral lines are well resolved. An ofi'set be- 



tween the image centre and the centre of the slit results in 
a change in the angle of incidence of light onto the grating, 
and thus a slightly blue or red shifted image. This causes a 
wavelength dependant change in position, and therefore an 
artificial spectroastrometric signature. As we use a wide slit 
of 5 arcsec this effect is likely to be origin of many of the 
artifacts present in our data. 

The empirical finding is that it is crucial to obtain mul- 
tiple spectra, comprising of anti-parallel sets of data. Such 
data will identify artifacts regardless of their cause, and can 
also be used to remove systematic efi'ects. 



4 SPLITTING BINARY SPECTRA 

As spectroastrometric signatures trace changes in fiux dis- 
tributions such signatures can be used to disentangle a con- 
volved binary spectrum into its constituent spectra. Two 
approaches can be used. One method relies upon a priori 
information while the other does not. 
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The first app roach was pioneered by iBailevI (|l998al ). 
and later used by iTakami et alj l|2003l ). The positional sig- 
nature of a binary system is directly proportional to the 
system separation and continuum flux ratio. Therefore, if 
these properties of a binary system are known, the intensity 
and positional spectra observed can be used to disentangle 
the individual fluxes of the two components. 

In contrast, the spectra splitting method of lPorter et al.l 
does n ot require any prior knowledge to separate bi- 
nary spectra. iPorter et al.l (|2004 ) present a series of simu- 
lations, in which the dependence of spectroastrometric ob- 
servables on the flux ratio and separation of a binary sys- 
tem are in vestigated. U s ing re lationships established by the 
models of lPorter et al.l (|2004l ). and the three spectroastro- 
metric observables (the centroid, total flux and width), the 
individual fluxes of the binary components can be recov- 
er ed. For the details of the method the reader is referred 
to lPorter et al.l (|2004l ). Essentially a model binary system is 
considered with a range of separations. For each separation 
the continuum flux ratio is estimated, from the observed 
width of the spectral profile, cr, in the continuum. Then, 
using the positional excursions observed, the binary o spec- 
trum is predicted. This is then compared to the observed a 
distribution. The best fit allows the binary separation to be 
estimated. Once the binary separation, and the associated 
continuum fiux ratio, have been deter mined the app roach 
used is essent i ally t he same as that of iBailevI (|l998al ) and 
iTakami et al.1 ll2003l). We discuss t he use of this method in 
more detail in IWheelwright et al . ('2009'). Here we attempt 
to apply the method of lPorter et al. (.2004) ) in the red region, 
as the Ha line is often associated with the largest features. 
We then use the det ermined proper ties of the binary system 
with the method of IBailevI jl998al ) to separate the binary 
spectra in the blue region 



4.1 Separated binary spectra 

To separate unresolved binary spectra into the constituent 
spectra it is required that prominent spectroastrometric sig- 
natures are observed across photospheric lines in the B re- 
gion. Also, we only attempt to separate component spec- 
tra when the spectroastrometric signatures trace a linear 
excursion in the XY plane, as opposed to a loop. If a spec- 
troastrometric signature is solely due to a binary system the 
signature will trace a linear excursion in the XY plane, as 
demonstrated by the example of GU CMa (Fig. [T|. There- 
fore, this criterion should exclude contaminated binary sig- 
natures and signatures not due to binary systems, issues 
discussed in Sections 15. II and 15. 21 As a result of these crite- 
ria it was not possible to separate the constituent spectra of 
all the binary systems detected. 

We separate the unresolved binary spectra of 9 systems 
into the constituent spectra (we present the binary prop- 
erties used/established in Table [3]). Spectral types for each 
spectra were determined by comparing the spectra to that 
of Morgan-Keenan standard stars and comparing ratios of 
key diagnostic lines. We present the results of assessing the 
spectral type of each component in Table |4] In addition, we 
determine the system mass ratios by assessing the mass of 
each component from its spectral type, using the data of 
iHarman^ (|l988h . In some cases the spectral types of each 



Table 3. A summary of the binary properties used/established 
when separating the unresolved spectra. Column 1 presents the 
binary systems for which the constituent spectra were separated, 
column 2 lists the binary separations used or established and col- 
umn 3 presents the binary separations in the literature. Column 
4 presents the binary PAs determined from the data discussed 
here, and column 5 lists the binary PAs from the literature. Fi- 
nally, column 6 contains the difference in brightness between the 
two components that was used. References for the literature val- 
ues are presented in Table [2] 



Binary 


d 




PA 




AB 




II 


II 


n 


n 


(magnitudes) 


HK Ori 


0.36 


0.35 


46.9 ±3.1 


41.8 


1.0 


T Ori 


0.84 




107.2 ±2.5 




2.5 


V586 Ori 


1.00 


0.99 


216.8 ±3.3 


30.3 


3.5 


HD 37357 


0.14 


0.19 


61.5 ±4.1 


49.0 


1.75 


V1788 Ori 


0.69 


0.52 


131.3 ±6.6 


352.9 


3.5 


HD 245906 


0.13 


0.13 


81.9 ±3.1 


77.1 


2.5 


GU CMa 


0.65 


0.65 


197.9 ±0.2 


194.5 


1.1 


MWC 166 


0.52 


0.65 


298.3 ±0.7 


297.8 


1.2 


11 Cep 


0.44 




54.3 ±2.0 




3.5 



component of a binary system had already been estimated, 
e.g. GU CMa and MWC 166. The spectral types deter- 
mined using the spectroastrometricall y split spectra are in 
good agreement with previous results (jBouvier fc CorporonI 
bOOlT ). This provides an important check on the validity 
of the spectroastrometric procedure. In addition, the spec- 
tral types determined for the primary components generally 
agree with previous classifications of the composite spectra, 
which also provides a consistency check. 

The separated spectra are presented in Fig. [5] From 
examination of the spectra (Fig. [2] and spectra split in the 
R band) it is clear that in some cases only the primary 
component is responsible for the emission lines seen in the 
comp osite spectrum. This is in agreement with the find- 
ing of lBouvier fc CorporonI (|200ll ). who report that in many 
HAe/Be binary systems only the primary exhibits a signifi- 
cant NIR excess, i.e possess circumstellar material. 

We compare the binary mass ratio distribution observed 
with that predicted assuming the secondary mass is drawn 
at random from the Initial Mass Function (IMF). For the 
determined mass of each primary we randomly draw a com- 
panion mass from the IMF given bv lKroupal (|200lh . To es- 
timate the most probable companion mass we do so 10,000 
times and use the resultant average mass. Table |4] compares 
the observed and the predicted mass ratio. The predicted 
mass ratio distribution peaks at a relatively low values, and 
no systems are predicted to have a mass ratio greater than 
0.1. In contrast, the observed mass ratio distribution is no- 
ticeably skewed towards higher values, see Fig. |3l 

We assess how different the two distributions are using 
the one sample Kolmogorov-Smirnov (KS) Test. Accord- 
ing to the KS test the scenario that the secondary mass 
is randomly selected from the IMF may be rejected with 
almost 100 per cent confidence. Thus it appears the mass 
ratio distribution of the binary systems is almost certainly 
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Table 4. The results of separating binary spectra into the two constituent spectra. Column 1 lists the objects in question and column 2 
denotes the spectral type of the primary. Column 3 contains the uncertainty in the spectral type of the primary while column 4 lists the 
spectral type of the secondary and column 5 presents the uncertainty in the spectral type of the secondary. The spectral types of these 
objects taken from the literature are listed in column 6. The resulting mass ratio is presented in column 7. Finally, the predicted mass 
ratio of the system, if the secondary were drawn at random from the IMF (see text for explanation), is listed in column 8. 



Binary 


Typei 


ATypei 


Types 


ATypea 


Spec Typeiit 


Qob 


^prcd 






(sub-types) 




(sub-types) 








HK Ori 


AO 


2 


K3 


3 


GlVeA,A4pevB 


0.33 


0.07 


T Ori 


A2 


1 


A2 


2 


A3IVevA,A0C 


1.00 


0.07 


V586 Ori 


A2 


1 


F5 


5 


A2V° 


0.65 


0.07 


HD 37357 


A2 


1 


A4 


2 


AOVe^ 


0.94 


0.07 


V1788 Ori 


A2 


1 


F5 


3 


B9Ve^ 


0.65 


0.07 


HD 245906 


Al 


1 


G5 


5 


B8eF 


0.52 


0.07 


GU CMa 


Bl 


1 


B2 


1 


B2vne<^ 


0.78 


0.02 


MWC 166 


BO 


1 


BO 


1 


BOIVeH 


1.00 


0.01 


11 Cep 


B3 


2 


B4 


2 


B2peH 


0.84 


0.03 



References: A) iMora et al.l 1I2OOII). B ) iBidelmar 
iHerbig fc Bell l|l988h . O lCuetted l|l968l '). HI iHiltnei 



\l954 ). C) [Hernandez et al.l ||2004|') . D) ISmithI lll972l') . E) iThe et al.l lll994f) , F) 
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Figure 3. The cumulative distribution of: the observed binary 
mass ratio {solid line) and the mass ratio distribut ion predicted 
by random sampling of the IMF, from lKroupal ll200lh {long dashed 
line). 



not determined by random sampling from the IMF. Clearly, 
this finding only retains its statistical significance if we can 
detect all mass ratios equally well. However, the lowest de- 
tectable mass is, on average, ~ O.OM©. This is as the sensi- 
tivity of spectroastrometry is limited by the relative bright- 
ness of binary components. If the primary component of a 
binary system is more than 5 magnitudes brighter than the 
secondary the system will probably not be detected by spec- 
troastrometry. As a result the lowest detectable mass ratio 
is ~ 0.3, greater than the location of the peak of the mass 
ratio distribution predicted. Therefore, it is possible that 
a large number of low mass ratio systems are undetected, 
which would introduce a bias to the mass ratio observed, 
skewing the distribution to high values. 

To quantify the effect this may have, we consider the 



case in which this bias has the largest effect possible. We 
assume every star that is a non detection is a binary sys- 
tem which we do not detect, due to a large difference in 
brightness between the two components. We assign each fic- 
tional system a mass ratio determined by the IMF, which is 
generally 0.1 and below our detection limits, and add these 
systems to our sample. Using the KS test we find that the hy- 
pothesis that the new 'observed' mass ratio originates from 
randomly sampling the secondary mass from the IMF may 
still be rejected with 99.55 per cent confidence. In this case 
the observed mass ratio distribution appears incompatible 
with the secondary mass being selected at random from the 
IMF at almost a 3a level. Therefore, the result appears ro- 
bust, despite our detection limit. 



5 DISCUSSION 

5.1 On the large FWHM features unaccompanied 
by positional features 

Many stars in the sample, such as AB Aur, present spec- 
troastrometric signatures in which the FWHM features 
are much more pro minent than any positional excursions. 
iBaines et al.l (|2006l ) suggest that these features are due to 
wide binary systems, where wide refers to a separation 
greater than half the slit width. In the case of AB Aur 
we detect a similar spectroastrometric signature over Ha 
to iBaines et al.l (120061 ). However, a 'wide' binary would be 
resolved in the data, as we use a slit of 5 arcsec. The longslit 
spectra were visually checked for evidence of a resolved com- 
panion, none was found. Therefore, as a resolved system does 
not create a spectroastrometric signature, there must be an- 
other, as yet unconsidered, source of the spectroastrometric 
signature. It is plausible that light from nebulosity could 
have distorted the spectroastrometric signatures. Extended 
emission is noticeable in many longslit spectra. However, it 
was found that masking the nebulosity had no effect on the 
spectroastrometric signature observed. 

As an alternate explanation we consider the sugges- 
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tion of iMonnier et al.l l|2006l ) that some spectroastromet- 
ric features could be cau sed by the pres ence of dusty ha- 
los around HAe/Be stars. iMonnier et all (2006) found that 
many HAe/Be stars, including AB Aur, are surround by ex- 
tended features of up to 0.5 arcsec. These features contribute 
up to 20 per cent of the NIR flux detected. Such halos are 
not well studied, but could constitute light scattered from 
the remnant natal envelopes of such stars, dust entrained 
in a wind or localised thermal em ission a few au from the 
central star l|Monnier et al.l [20061 ). Such extended emission 
would be unresolved in a longslit spectrum, and could lead 
to an increase in the FWHM while not changing the photo- 
centre position. However, this requires that the line profile 



of the scattered light is differe nt from the original e mission 
source profile. As discussed by IMonnier et all (l2006l ) this IS 
certainly plausible. A non-uniform distribution of Ha fiux 
and line of sight dependent absorption would both result in 
the observer and the scattering media seeing slightly differ- 
ent line profiles. 

We explore whether an unresolved, extended halo could 
result in a spectroastrometric signature, similar to that ob- 
served over the Ha line in the case of AB Aur and other 
stars, using a simple model. The model treats a single star as 
a point source and surrounds the star with a halo which con- 
tributes 20 per cent of the total flux. Here the halo is offset 
from the central star position by approximately 0.50 arcsec. 
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Figure 4. The input spectra, the resultant spectrum and the 
spectroastrometric signature of the model comprising a star sur- 
rounded by a halo. The halo is modelled as a uniform ring from 
0.45 arcsec to 0.60 arcsec, centred on a point 0.075 arcsec offset 
from the star. The seeing used was 1 arcsec and the halo con- 
tributes 20 per cent of the total flux. The long dashed lines are 
the spectroastrometric signature of AB Aur over the Ha line, in 
the North-South direction and the short dashed lines are the spec- 
troastrometric signature of AB Aur in the East- West direction. 
The dot- dashed line is the averaged spectrum of AB Aur in the 
Ha region (the AB Aur data has been rescaled in the dispersion 
direction). 



The flux emanating from the halo has a uniform distribu- 
tion in space. A P-Cygni type profile is assigned to the star 
and a similar line profile, minus the blue-shifted absorption 
component and with a slightly different line to continuum 
ratio, assigned to the halo flux. The total flux distribution is 
mapped onto an array representing a CCD chip. The array 
is then convolved with a Gaussian in the spatial direction to 
represent the effects of seeing. Finally, the output spectrum 
is extracted in a standard fashion and spectroastrometry is 
conducted on the artificial observation. We present the re- 
sults of this exercise in Fig. |3] 



There exists a qualitative similarity between the model 
and observed signatures (Fig. |4} . The model did not com- 
pletely recreate the extent of the FWHM feature observed in 
the case of AB Aur. However, given the unknowns involved, 
e.g. the amount of light scattered and the extent of the halo, 
this does not exclude this scenario. Therefore, we conclude 
that it is likely that FWHM features accompanied by small 
or nonexistent positional signatures are due, at least in part 
to unresolved, extended, halos. Alternatively, a wind could 
also res ult in a similar positio nal spectroastrometric signa- 
ture, see lAzevedo et al] l|2007l ). This has important implica- 
tions on splitting the binary spectra, which are discussed in 
Section 15.21 We note that some known binary systems ex- 
hibited larger FWHM features than positional features, and 
as such this is not a unique diagnostic. 



5.2 An ev aluation of the method of IPorter et all 
(|2004l l 

Implicit in both methods of splitting spectra is the assump- 
tion that the system in question comprises of two point 
sources. In Section 15.11 we demonstrat e that dusty halos , 
which surround som e HAe/Be stars l|Leinert et al.l 1200 ll : 
iMonnier et al.|[2006l ). can give rise to spectroastrometric sig- 
natures. It may be expected that if the spectroastrometric 
signature of a binary system is contaminated by the signa- 
ture o f a halo, the spectra splitting method of lPorter et al.l 
(|2004 ) will not be able to correctly separate th e constituent 
spectra. In many situations the method of IPorter et al.l 
( 2004) failed to fit the observed FWHM spectrum of a known 
binary system, even when the separation considered was in- 
creased to many times the binary separation. Here we in- 
vestigate whether this could be due to contamination of the 
binary spectroastrometric signature by an additional, unre- 
solved source of fiux. 

We construct a model of a binary system with a sepa- 
ration of 0.3 arcsec. The binary system has a difference in 
brightness of 3.5 magnitudes and is surrounded by a halo 
that extends from 0.4 arcsec to 0.6 arcsec from the central 
star. An artificial longslit spectra is generated and spectroas- 
trometry is applied to the synthetic data to generate the ob- 
servables necessary to separate the cons tituent spectra. Fi- 
nally, the method of lPorter et al.l (|2004l ) is used to attempt 
to split the unresolved binary spectrum into its constituent 
spectra. 

The results of first modelling the aforementioned binary 
system without an extended halo component, and the results 
of including an extended halo component, are displayed in 
Fig.O When no halo component is added the two spectra are 
clearly separat ed, demonstrating t he power of this approach. 
The method of lPorter etHI (|2004l ') fits the observed FWHM 
spectrum, and as a consequence splits the binary spectra 
correctly. In contrast, when the h alo component is a dded 
to the binary model, the method of lPorter et al.l (|2004| ) can 
no longer correctly sep a rate t he constituent spectra. The 
method of IPorter et al.l (|2004l ) no longer fits the observed 
FWHM signature, as shown, and consequently fails to sep- 
arate the two binary spectra correctly. This is only to be 
expected as: a) the positional and FWHM features observed 
are no longer due to two point sources and b) the method 
attempts to apportion the observed fiux, which is due to 
three sources, to only two sources. 

This would also be the case if the spectroastrometric 
signature observed were due to a triple system. The degree 
to which a third component would compromise the spectra 
splitting procedure would depend on the relative brightness 
of the system components. For example, the least bright 
component would have to be brighter than 1 per cent of 
the combined fiux emanating from the two brightest com- 
ponents to contaminate the spectroastrometric signature. A 
triple system might be expected to exhibit distinctly differ- 
ent spectroastrometric signatures over different lines. The 
norm for this sample is for the spectroastrometric signatures 
over different lines to be consistent, as demonstrated by the 
example of GU CMa (Fig. [!}. Therefore, if triple systems 
are present in the sample it would appear that the tertiary 
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Figure 5. The results of splitting the spectra of the model of a binary system (leftmost two panels) and the model of the same binary 
system plus an unresolved halo which contributes 20 per cent of the system flux (rightmost two panels). In the panels presenting the 
o bserved a dist r ibutio n the solid line is the observed quantity while the dashed line is the best fit a distribution predicted by the method 
of lPorter et"al] ll2004h . 



components are not bright enough to significantly effect the 
spectroastrometric signatures observed. 

In summary, the method of lPorter et al. I (|2004l ) is com- 
promised if an additional source of flux, besides the binary 
system, is present. Such a source may a tertiary stellar com- 
ponent, a dusty halo, or material in a wind. We suggest that 
this is the reason t hat, more times than not, the method of 
IPorter et all (|2004l ') clearly does not fit the observed FWHM 
features and, as a consequence, fails to separate the spectra 
of many unresolved binary systems. If this is the case for 
the systems where we apply the spectra splitting procedure, 
the returned spectra will not be correctly separated. How- 
ever, we only present separated spectra for systems whose 
spectroastrometric signature appears solely due to a binary 
system, with none of the complications mentioned above (see 
Section [4. II for the condition used to exclude contaminated 
signatures) . 

5.3 On the separation of HAe/Be binary systems 

The previously detected systems in the sample have physical 
separations between ~ 40 and 1200 au. Unfortunately, this 
cannot be easily translated into a separation distribution as 
the stars are at very different distances (between 143 and 
2000 pc), and also due to the various selection effects in dif- 
ferent detection methods. Of the newly discovered binaries 
the angular separation to which we are sensitive is in the 
range 0.1 - 2.0 arcsec. However, the different distances to 
each star change the physical separation to which this cor- 
responds. For the nearest system at 143 pc, the separation 
range is 14 - 285 au. For a system at the average distance 
of a star at 600 pc it is ~60 - 1000 au. For the most distant 
system at 2000 pc, it is 200 - 4000 au. 

None of the previously detected binaries are closer than 



30 au, which is the peak of the field G-dwarf distribution. 
Half o f all G-dwarf binary systems have separations less than 
30 au (|Duauennov fc Mavor|[l99ll ) . Of the newly discovered 
binaries, only one could be closer than 30 au (V1366 Gas 
at 164 pc), and two other systems could have separations as 
small as 30 - 40 au. Therefore, at least 60 per cent of HAe/Be 
stars have a companion between about 30 - 4000 au, and 
probably in the range 60 - 1000 au. This is significantly 
greater than the fraction of G-dwarfs at the same separa- 
tions (around 40 per cent between 30 - 4000 au, and 25 
per cent between 60 - 1000 au). This overabundance of bi- 
naries is not dissimilar to, but apparently larger than, the 
overabun dance of binary syste ms found in young T Tauri 
stars fsee iDuchene et all (|2007l ) and references therein). 

Thus, unless there is an almost complete lack of com- 
panions < 30 au, it is difficult to imagine that the binary 
fraction of HAe/Be stars is much less than 100 per cent. If 
HAe/Be stars exhibit a similar abundance of companions 
< 30 au to G-dwarfs this would suggest that many HAe/Be 
stars are triple or higher-order systems. Indeed we present 
the detection of four additional components in previously 
detected binary systems, meaning that these systems are at 
least triple systems. 

It is worth noting that many of these systems are rela- 
tively soft, with separations greater than a few hundred au. 
As a result, these sys tems are susceptibl e to destruction 
in dense clusters (see IParker et al.l (|2009l ) and references 
therein). This suggests that many of these HAe/Be stars 
have not spent a significant time in a very dense environ- 
ments, e.g. densities of > 10^ M0 pc~^, which are not 
unusual in star fo rming regions of young clusters. Indeed, 
iTesti et all l|l999D found that no Herbig Ae stars are asso- 
ciated with clustered environments, and that while Herbig 
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Be stars are sometimes situated in a small cluster, the asso- 
ciated stellar densities are approximately 10^ — lO^pc"^. 

5.4 The mass ratio and formation mechanisms of 
H Ae/Be binary systems 

It appears that the mass ratio of Herbig Ae/Be binary sys- 
tems is skewed towards relatively high values, and is incon- 
sistent with random sampling from the IMF. However, as 
the spectra splitting technique did not work in every case, 
the sample size is too small to attempt to constrain the un- 
derlying distribution of mass ratios. Instead, we discuss the 
possible implications this finding has on the formation mech- 
anisms of intermediate mass stars. We note that random 
pairing of binary components has already been exclude d 
in the OB association Sco OB2 (|Kouwenhoven et allliooH ). 
Here we extend this finding to younger systems, and higher 
masses. 

iBaines et al.l l|2006l ) found that the circumstellar discs 
of the components of Herbig Ae/Be systems are preferen- 
tially aligned with the binary position angle. This already 
sug gests that the second ary formed by disc fragmentation. 



: Goodwin et al.l |20o3). As noted bv iKouwenhoven et al.l 
l|2009t ). disc fragmentation would be expected to produce 
stars of roughly similar mass (within a factor of a few). 
Disc fragmentation should occur during the earliest phases 
of star formation. During such phases there is an abundance 
of gas to accrete, and the circumprimary disc is still mas- 
sive enough to fragment. The secondary in the disc is able 
to accrete material from the disk more easily than the pri- 
mary, as the angular momentum o f the material is closer 
to the secondary than the primary (jWhitworth et all 1 19951 : 
iBate fc Bonnelllll99'iD . Therefore, this scenario results in a 
binary system with a high mass ratio, higher than if ran- 
dom sampling from the IMF determined the mass of the 
secondary. Indeed, recent models of massive star formation 
demonstrate that binary systems with high mass ratios (0.7) 
and large separ ations (^ 1000 au) ca n be formed from disk 
fragmentation (iKrumholz et al.ll2009t ). 

The separations of the binary systems in our sample also 
suggest disc fragmentation as the mode of binary formation. 
IWhitworth fc Stamatello^ (|2006l ) show that a massive disc 
can fragment beyond a critical radius -Rfrag which depends 
on the mass of the primary Mi, as follows: 

M ^ 
R... > 150 (^^ 

For a Herbig Ae/Be star ~ lOM© , and so 7?frag > 300 au, 
which is a typical separation of the systems in our sample. 

Therefore, the properties of the Herbig Ae/Be binary 
systems observed indicate that these systems formed via 
disk fragmentation. Given that the sample includes stars as 
massive as ^ 15M0, this favours the core collapse and sub- 
seque nt monolithic accretio n scenario of massive star forma- 
tion (IKrumholz et al. 20091). as opposed to the merger and 



captu re scenarios ( Ballv fc Zinneckerl2005l : lMoeckel fc Ballvl 
I2OO7I ). 

The Herbig A e/Be stars in this sample are not located 
in dense clusters (jTesti et al.lll999l ). In addition, the wide 
separations of the binary systems and their young ages sug- 
gest that they formed in isolation. Firstly, binaries this wide 



are relatively soft and could not hav e spent a significan t 
amount of time in a dense cluster, see IParker et all (|200i). 
Secondly, no binary this wide could have survived ejection 
from a cluster. Together this suggests that a fairly large core 
(> IOM0) formed in relative isolation and produced a mas- 
sive binary system, rather than a small cluster. That these 
H Ae/Be stars formed in a massive, isolated core shows that 
competitive accretion, e.g. iBonnell et al.l (| 19981 ). is not re- 
quired to form stars of up to at least IOM0, as presumably 
no larger reservoir of gas existed beyond the single core. 

This cannot be infrequent as we find several Herbig 
Ae/Be systems which fit this pattern. However, the sam- 
ple we have used is not complete, the population of HAe/Be 
is heterogeneous to begin with and the selection criteria may 
well impose certain selection effects on membership of the 
HAe/Be class. Biases and incompleteness are impossible to 
fully quantify, but we can state that a not insignificant frac- 
tion of A/B stars can form in isolation from a massi v e core . 
This is in qualitative agreement with Ide Wit et all (|2005h . 
who repo rt that even O type stars may form in isolation. In 
addition, IParker fc Goodwin! (|2007l ) also find that a few per 
cent of massive stars might form in relative isolation. 



6 CONCLUSIONS 

In this paper we present spectroastrometric observations of 
a relatively large sample of HAe/Be stars. Here we present 
the salient findings of this work: 

• We find a high binary fraction, 74 ± 6 percent, consis- 
tent with previous studies. 

• Using spectroastrometry to separate the unresolved bi- 
nary spectra we determine spectral types for the components 
of 9 systems. 

• The mass ratios of these systems, determined from the 
constituent spectral types, are inconsistent with a secondary 
mass randomly selected from the IMF. 

• Although our sample is small this result constrains the 
mode of binary formation in that the mass ratios and sep- 
arations of the binary systems observed suggest that the 
secondary forms via disk fragmentation. 

• The properties of the binary systems observed indicate 
that these systems have not spent a significant amount of 
time in dense, clustered environments. Therefore, these sys- 
tems demonstrate that isolated star formation can produce 
stars as massive as ~ 10 — ISMq. 
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APPENDIX B: HI SPECTROASTROMETRIC 
SIGNATURES 

In Fig. IBll the Hq or H/3 profiles and the associated spec- 
troastrometric signatures of the 47 stars in the sample are 
presented. For each object the average, normalised, intensity 
spectrum is presented alongside the position spectra in the 
two perpendicular directions observed. In addition, we also 
present the FWHM spectra, also in these two directions. To 
keep the appendix concise only one profile per observation is 
presented, where possible the Hq signature, and where not 
the II/9 signature. Artifacts are indicated by dashed lines. 



APPENDIX A: A SUMMARY OF THE 
SPECTROASTROMETRIC SIGNATURES 

In Table IXn we present a summary of the spectroastromet- 
ric signatures over the Ha or H/3 lines for the 47 stars in 
the sample. Where possible the signature over Hq is sum- 
marised. However, not all the objects in the sample were 
observed in the R band. Therefore, in these cases the prop- 
erties of the H/3 signature are presented. 
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Table Al. A summary of the spectroastrometric results across H« or H/3. Column 1 lists the objects observed, columns 2 and 3 contain 
the continuum uncertainty in the position and FWHM spectra respectively, column 4 lists the average equivalent width of the line in 
question (accurate to 10% on average), column 5 denotes the emission profile type while columns 6 and 7 list the observed change in 
centroid position and FWHM over the line, and column 8 contains the calculated PA of the systems detected. 
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FWHM a 
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Ha profile^ 
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AFWHM 


PA 
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(A) 
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(o) 
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97 ± 3 
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3.0 


6.8 
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IIR 


Artifact 


61 ± 7 




IP Per 


1.8 


4.4 


2.8 


M 


<8 


15 ± 5 




AB Aur 


1.2 


2.6 


—22 


IVB 


14 ± 1 


233 ± 3 


45.8 it 5.1 


MWC 480 


1.1 


2.9 


— 18 


IVB 


^4 


145 ± 3 




UX Ori 
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3.0 
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HD 35929 
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t Profile classification from iReipurth et al.l ([199^ (I: symmetric emission, II: double peaked emission where the secondary peak is 
greater than half the intensity of the primary peak, HI: double peaked emission where the weaker peak is less than half the intensity 
of the stronger peak, IV: P-Cygni profile. The position of the weaker peak (or absorption component) with respect to the central 
wavelength is indicated by R or B. Absorption profiles are designated by Ab. Profiles with multiple absorption components are 
designated by M.) 
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Figure Bl. Ho profiles and spectroastrometric signatures. From left to right: VX Cas (data from the WHT), VX Cas (data from tiie 
INT), IP Per, AB Aur, MWC 480, UX Ori, HD 35929, V380 Ori and MWC 708. 
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Figure Bl. H« profiles and spectroastrometric signatures. From left to right: HD 244604, HD 37357, MWC 120 (data from the WHT), 
MWC 120 (data from the INT), MWC 137, HD 45677, LkHa 215, MWC 147 (data from the WHT) and R Mon. 
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Figure Bl. Ho profiles and spectroastrometric signatures. From left to right: V742 Mon, GU CMa (data from the WHT), MWC 166, 
HD 76868, MWC 297, HD 179218, HD 190073, BD+40 4124 and MWC 361 (data obtained with a 2.5 arcsec sht). 
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Figure Bl. Ha and H/3 profiles and spectroastrometric signatures. From left to right: MWC 361 (data obtained with a 4 arcsec slit), 
II Cep, BHJ 71 (data from the WHT), MWC 1080, V594 Cas, V1185 Tau, V1012 Ori, V1366 Ori and V346 Ori. 
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Figure Bl. H/3 profiles and spectroastrometric signatures. From left to right: HK Ori, V1271 Ori, T Ori, V586 Ori, V1788 Ori, HD 
245906, RR Tau, V350 Ori and MWC 790. 
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Figure Bl. H/3 profiles and spectroastrometric signatures. From left to right: MWC 147 (data from tlie INT), V590 Mon, OY Gem, 
GU GMa (data from tlie INT), HD 81357, SV Cep, MWC 655 and BHJ 71 (data from the INT). 



